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Abstract 

In this paper we show how polynomial mappings of degree K from a union of disjoint intervals onto 
[—1,1] generate a countable number of special cases of generalizations of Chebyshev polynomials. We 
also derive a new expression for these generalized Chebyshev polynomials for any genus g, from which 
the coefficients of x" can be found explicitly in terms of the branch points and the recurrence coefficients. 
We find that this representation is useful for specializing to polynomial mapping cases for small K where 
we will have explicit expressions for the recurrence coefficients in terms of the branch points. We study in 
detail certain special cases of the polynomials for small degree mappings and prove a theorem concerning 
the location of the zeroes of the polynomials. We also derive an explicit expression for the discriminant 
for the genus 1 case of our Chebyshev polynomials that is valid for any configuration of the branch point. 
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1 Introduction and Preliminaries. 

Akhiezer PP and, Akhiczer and Tomcuk [3j introduced orthogonal polynomials on two intervals which 
generalize the Chebyshev polynomials. He observed that the study of properties of these polynomials requires 
the use of elliptic functions. In the case of more than two intervals, Tomcuk |17|. investigated their Bernstein- 
Szego asymptotics, with the theory of Hyperelliptic integrals, and found expressions in terms of a certain 
Abelian integral of the third kind. However, in his formulation, certain unknown points on a Hyperelliptic 
Riemann surface emerge due to the lack of an explicit representation of the original polynomials. This was 
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circumvented in |^. In his book on elhptic functions 4 Akhiezer, obtained exphcit formulas for the two 
interval case as an example of the application of elliptic functions to the theory of conformal mapping. 

In 1984 Al-Salam, AUaway, and Askey introduced sieved ultraspherical polynomials which arc orthogonal 
with respect to an absolutely continuous measure supported on [— 1 , 1] but the weight function vanishes at 
fc + 1 points. Ismail jl2| observed that the vanishing of the weight function means that the polynomials 
are orthogonal on several adjacent intervals. He then introduced one additional parameter in the definition 
of the sieved ultraspherical polynomials which made them orthogonal on several intervals with gaps. In 
particular his polynomials include analogues of the Chebyshev polynomials of the first and second kind. 
Their continuous spectrum is TkW-c.c]), for c e (-1, 1), where is a Chebyshev polynomial of the first 
kind. More over these polynomials have simple closed form expressions and elementary generating functions. 
More general sieved polynomials are in f^. 

In jl4| a study of Chebyshev type polynomials was carried out. They arose as extremal polynomials and 
are orthogonal on several intervals. The Chebyshev type polynomials satisfy similar extremal properties to 
the classical Chebyshev polynomials on [—1,1]. In particular they deviate least from zero on a union of 
disjoint intervals. The extremal polynomials also have the property that they are orthogonal with respect 
to some weight function of the form, 



rg-i 



[x] 



X e E, 



where E :— U^f^^^(aj, aj+i) and the polynomial 7rg_i(x) is chosen such that the weight is positive on the 
interior of E. The existence of an extremal polynomials on any set E is clearly not guaranteed. It was shown 
in |14| that the existence of an extremal polynomial on a set E is equivalent to the asymptotic periodicity 
of the recurrence coefficients of a sequence of orthogonal polynomials on E. The extremal polynomial also 
has the property that it maps the set E onto [—1, 1]. 

In |ll)j. using a polynomial mapping of this kind, a prescription was given for the construction of a sequence 
of orthogonal polynomials on E from a sequence orthogonal on a single interval. The corresponding weight 
supported on E can be obtained using the following procedure, outlined in |1U| . In this generality closed 
form expressions and generating functions are not available. This raises the question of finding a class of 
polynomials as general as possible which has the extremal properties of the polynomials in Peherstorfer 
and the simple structural formulas of Ismail's polynomials in jl2| . The first step was carried out in 9 and is 
further developed here. The present paper we develop this work further by borrowing ideas from polynomial 
mappings and sieved polynomials. 

We start by choosing a polynomial of degree ^ — 1 and call this polynomial W^_i (i). Then we form the 
expression. 



(II) ^^fj-iy-) ^ yyA-iyi-) ^ 

where, the M~^{x) are a complete set of inverse branches of Mj^{x) and Ifi is the leading coefficient of Mf^(t). 
Rearranging the above equation and taking the limit as i — + M^^{x) we find, 

^ W,.,iM-\x))^ 
M'^{M-\x)) 
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By comparing coefficients of t in (1.1) it is clear that, 
J? 

^UJj{x) = c 

where c is the quotient of the leading coefficients of T4^^_i(a;) and M^{x). So, having chosen a particular 
polynomial Wii-i{x) we are in a position to define the transformed weight function in the following way. 
Let Eq be the interval [a, 6], and Pn{x) the polynomials orthogonal with respect to wo{x) over Eq : 

hnSnm = / Pn{x)P„i{x)wo{x)dx 
J a 

= Pn{x)Pm{x) ^^LJj(x)| WQ{x)dx 

-E / Pn{M^{t))P^{Ms,{t))^^^^,^M'^{^^^^^ 



c 

1 

c Je 



Pn{MsK{t)) Prn (M^ it))WA^iit)wo{t)dt, 



where 



and, 



E := uf^.Aj, 



: [a,b]^Aj. 

So the new sequence of polynomials is orthogonal on the disjoint set E with respect to the weight function 

w{t) = Wsi-i{t)wo{Msi{t)). 
It was shown in ^U] that 

Pndx) = PniM^ix)) n = 0,l,2,... 

This generates every polynomial in the sequence. The polynomials Pn^+j{x) I < j < A — I, are 
referred to as the intermediate polynomials. Formulas relating the intermediate polynomials and recurrence 
coefficients of the new sequence to known quantities of the old sequence are given in . 

Some explicit expressions for the Chebyshev polynomials on several intervals, along with many other 
important results can be found in 0. This paper is organized as follows. In this section we recall the theta 
function representation of these polynomials which was derived in We also present some facts concerning 
the polynomials that will be referred to throughout this paper. 

Section 2 is a study of two auxiliary polynomials Sg{x;n) and Gg+i{x;n), of degree g and 5 + 1 re- 
spectively. The auxiliary polynomials are fundamental to the theory of the generalization of Chebyshev 
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polynomials under consideration. We also derive for the first time non-linear difference equations satisfied 
by the recurrence coefficients that appear in the three term recurrence relation. 

In section 3 we derive a product representation for our polynomials from which one can obtain the coef- 
ficients of the polynomial in terms of the recurrence coefficients and the branch points. This representation 
is a new one that involves the auxiliary polynomials studied in section 2. It is particularly useful when 
we consider certain configurations of the branch points where the recurrence coefficients simplify into an 
algebraic form. 

In section 4 of this paper we show how certain polynomial mappings give rise to a countable number 
of special cases of our generalization of the Chebyshev polynomials. These cases correspond to when the 
recurrence coefficients are periodic. We give explicit formulas for the polynomials of degree not equal to an 
integer multiple of the period. 

The monic form of the Chebyshev polynomials on several intervals will be denoted by They are 

orthogonal on E := [—1, a] U^Zj [Pj, aj+i] U [(3g, 1] with respect to the following weight, 



(1.2) 




so that, 




The polynomials of the second kind are 



(1.3) 




Both Pn{x) and Qn{x) satisfy the recurrence relation 



(1.4) 



(x) = Un+i{x) + bn+iUn{x) + a„M„_i(a;), 



see 12]. The Stieltjes transform of the weight function is 




Note that. 



ip{x) 



—w{x). 



In Chen and Lawrence considered a hyperelliptic Riemann surface of genus g, 



a 



z-aj){z- (3j), 
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from which a representation for the polynomials in terms of the associated Riemann theta function was 
constructed. They also gave a representation for the recurrence coefficients in terms of the theta functions. 
To formulate these representations the following facts regarding this Riemann surface are required. 

A canonical basis of cycles is chosen on the surface that is composed of 2g closed contours. We denote by 
the a-cycles the g closed contours that lie on one sheet of the surface, and the b-cycles to be the remaining 
g contours that travel over both sheets. For a surface of genus g a basis of holomorphic differentials can be 
written in terms of g linearly independent differentials, 

dx X dx x^^^dx 



y{x)' y{x)'""' y{x) 
If we let, 

then the dujj defined as, 

form an orthonormal basis such that, 
duJk = 6jk- 

We can then define the period matrix B whose entries are, 
Bjk — I dujk- 

The Riemann theta function associated with the surface is defined to be, 
■d{s; B) ^ z9(s) ^ exp{iTr{t, Bt) + 2ni(t, s)), 

where s is a g-vector, and (m, v) denotes the standard scalar product of two vectors. Convergence of the 
series is assured since 3-B is a positive definite matrix. The theta function has the following properties. For 

t e Zf, 

di-s) = z9(s) 

Both the polynomials and the recurrence coefficients are given below. The expressions are taken from [H]. 
However, in view of lemma 10.2 in [S] and the quasi-periodicity of the theta functions we have simplified the 
expressions slightly. 



(1.5) Pn{x) = ^ (f«(p.)+4(px)) 
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where, for n > 1, 



and, 



^^(p^) = 2^-nto+n(p.))^(u+)^(u +nB) 



^ ^ i?(u+ - nB)^(u^) 



Using this notation. 



(1.6) g„(x) = ^ - f„(Px)) 



Here, pa; is a point on the top sheet of the Riemann surface that corresponds to a point x in the complex 
plane and, 



- 1 



pi 

dt. 



y{t) 



The kj can be determined from the condition that, 
I dfl = 0, j = l,...,g. 

Jaj 

Uj; is defined to be a g- vector, 

U:, := i doJi, / dojg] . 

and u"*" := Uoo± . If B is the vector whose components are, 
1 



Bj = — / d^, j = l,...,g. 



then, 

B = u+ - u" (mod Z9). 
Finally, the equilibrium potential of the set E is 



.00+ / 1 



dt. 



y{t) 

With this notation, we have the following expressions for the recurrence coefficients, 
2e-2xo]!i4±|i ifn = l 

^ ' " » -2Yn »?(u+-(n-l)B)^(u++»B) 

i?(u+ + (n-l)B)iJ(u+-nB) 
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and, 



(1.9) 



Here, 



i9(u+ -B) i9(u+) i9(u+-(n-l)B) i>(u+-nB) 



Also, 



(1.10) 



/.„ = 2e--"xo^(u" + »B)_ 
i?(u+ - nB) 



Theorem 1.1 T/ie polynomials of the second kind Qn{x) satisfy the following orthogonality relation, 



(1.11) 
where, 



?«(a;)Qm(x)— -- 



dx — ^nmhnt 



hn = n^hn- 



Proof : We wish to evaluate the fohowing integral, 

Qn{t)^ 



W{t) 



■t^dt 0<j<n-l. 



We have. 



Qn{z) = Pn{z)Hz) 



w{t)dt 



= Pn{z)i^{z) --J 

z Je^- 

= p^{z)^{z)-- [ p^{t)f2(-ydt 

= Pn{zMz)- 



where the last line fohows from the orthogonahty of the P„ (x) . Therefore, 



Quiz) 



Pn{z) 



1 



'4){z) \z 



hr, „ / 1 

+ 



yn+1 



n+2 
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We rewrite Ij as a contour integral using a closed contour A 9 [Fig. 1], in the slit complex plane that 
encircles all of the branch points. By continuous deformation of the contour onto the intervals that make up 
E we have, 



' Je w{t) 

So Ij can be written as, 

Ja Ja Hz) 

The differential Pn{z)z^dz has no residue at infinity for any j. The differential ^n+^'^f^-j dz has no residue 
at infinity for j — 0, ...,n — 2. Consequently, for these values of j, both integrals are zero which proves the 
orthogonality of the Qn{z) with respect to The integral /„_i is equal to the square of the weighted L2 

norm of Qniz)- In this case there is a residue at infinity. The residue is —/in and therefore, 

□ 

Theorem 1.2 Let E be the reflection of E around the point a; = 0, namely 

E := [-1, -fig] U]^g_i U [-ai, 1]. 

Let {Pn(x)} he the monic polynomials orthogonal on E with respect to w{x), where, 

, 1 / uu^^+Pj) 

^^^^ {l-x^)UU(^ + ^^y 

Then we have for n> g 

Pn+giPl) 

Pn+g{l3g) 

Qn+g{0Ll) ■ 

Qn+g{ag) 

Pn{-X) ... Pn+gi-~x) 

If n < g then, 



Pn-g{0l) 



(1.12) P^ix) 



(-l)"+f 



Pn-g{f3g) 
Qri-g(ai) 



Qn-g{0ig) 
Pn^g{~x) 
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Pnix) 



(-1) 



n+g 



Po(/5l) ... Pn+giPl) 



PoiPg) 



a-! a\ 



a? " ^ '9o(ai) 



Pq{-x) 



Pn+g{Pg) 
Qn+g(ai) 

Qn+g{ag) 
Pn+gi-x) 



In both equations, Dn is equal to the coefficient of Pn+g{—x) appearing in each determinant. 
Proof: Starting with the orthogonahty relation of the Pn{x) we make the change of variables x —x. 



_ [ ~ ~ I nj=i(2;-«j) 

= P^i-x)Pn{-x)^ ^^^^ 



So the Pn{—x) are orthogonal with respect to a weight function that is w{x) multiplied by a rational 
function. In general if we know the polynomials orthogonal with respect to a weight and we wish to 
find the polynomials orthogonal with respect to R{x)fi{x) where is a rational function that is positive 
on the interval of orthogonality we may invoke a theorem of Uvarov that gives these polynomials in terms 
of the original ones. The theorem was first published in jTJl and also appears in In our case the 

rational function has zeros and poles at the end points of the subintervals that make up E. However, a 
brief examination of the proof will show that the theorem is still applicable in this case. The result follows 
directly from the theorem. □ 

In section 2 we wish to study the following polynomials, 

(1.13) Sg{x;n) f[{x - a,) - Ql{x){x^ - 1) f[{x - (3,) 
and, 

a 9 

(1.14) Gg+i{x; n) P„(x)F„_i(a;) Y[{x - a,) - Qn{x)Qn^i{x){x^ - l)Y[{x - /3,). 

These functions are fundamental in the theory of the generalizations of Chebyshev polynomials considered 
here. 
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Theorem 1.3 The functions Sg{x; n) and Gg+i{x; n) are polynomials of degree g and g + 1 respectively. 
Proof: Prom (1.3) we know that, 

Qn{z) = Pn{zmz)- [ ^w{t)dt. 

Expanding the integrand in powers of t we have, 
1 1 i ^2 

7 = - + — + — + •••• 

z — t Z Z'^ 

Therefore, as ^ ^ oo+, we have 

P„{z)i;iz) - Qn{z) = O 
and similarly, 

Pn{z)ij{z) + Qn{z) = 0{z^-^). 

Multiplying the two together gives 
P^{z)il;\z) - Qliz) = O 
Therefore, 

g g 

l[{z - aj)P^{z) - iz' - 1) - mli^) = O (z^) . 

Since the left hand side is a polynomial it follows that it must be of degree g. For (1.14) we form the following 
product, 

(P„(Z)V(^) - Quiz)) {Pn-l{z)ij{z) + Qn-l{z)) 

= Pn{z)Pn-liz)'4){z) - Q„{z)Qn-l{z) - tp{z)hn-l 

where we have used, 

Pn-l{z)Qn{z) - Pn{z)Qn-l{z) = K-l, 

which is the Wronskian. As a result, 

g g 
JI{Z - aj)Pn{z)Pn-l{z) - {Z^ - 1) ^[(^ " Pj)Qn{z)Qn-l{z) = O {z^+^) . 

□ 
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From (1.13) and (1-14) we have 

(1.15) £n{px)£n{Vx) - 

Also, 



Sg{x]n) 
YVj=i{x-aj)' 



(1.16) £„(p.)£V,(p.) = Gg,.{.;n)^K^M^) 



(1.17) 5n(pa;)f»-l(px) 



Gg+i(x\n) ~ hn^iy{x) 



n?=i(^"".) ■ 

An alternative representation for the polynomials of the first and second kind was also given in |^ 
resulting from a coupled system of differential equations. 

(1.18) P„(z) := p„(^)cos1'„(z) 
and, 

(1.19) Qn{z) := iiP{z)pn{z) sin-^niz), 
where, 



Pn{z) = 



\3 



and, 



Ji y 

where the branch of the square root is chosen in such a way that Pn{z) — > V27i^, as 5Rz oo and 2/i„ is 
the leading coefficient of Sg{z\ n), to be shown later. Here the path of integration is from 1 to an arbitrary 
complex point z with the property 3z > 0, and the path of integration is entirely in the upper half plane. 
Expressions for the real values are then obtained by the analytical continuation of those above allowing 

2 ^ X e R. 

In these expressions, the 7j(n) are the zeroes of Sg{x, n), and Ui = y{px)\f^=-yj{n)- The Cj{n) are the first 
5 + 1 expansion coefficients of ^ \n£n{px) in terms of the following basis, 

1 a; a;'' 1 1 U + TJi U + Ug 



y' y''"' y ' X - ai' '"' X - ag' y{x - Jiin))''"' y{x - ^g{n)) 

This is the basis for meromorphic functions on the Riemann surface defined by y{x) whose zeroes must 
include simple zeroes at ±oo and whose poles must be from amongst the set pi, p_i, pQ,^ , pp. , p^^ : j = 1, .., g, 
where the poles at the a points can be at most double poles and the other poles must be simple poles. 
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2 Evaluating the polynomials Sg{x;n) and Gg+i{x; 



n 



In this section we give an algorithm for determining the coefficients of Sg{x;n) and Gg+i{x;n) in terms of 
the recurrence coefficients and the branch points. Let 



g 

J 



and, 



Sg{x;n) = J2vj{n)x- 

3=0 



9+1 

Gg+i{x;n) ='^^j{n)x^. 



Then we have, 



(2.1) mP^iz) - "^"^'^ - + - + 



ip{z) y{z) 



(2.2) ^(^)P„(^)P„_i(^) - -^^ , 

To find the r?j(n) coefficients we integrate both sides of (2.1) along the contour A. Note that, 
i,{z) = '^W{z), 

where 



W{z) = -. 



TV(^'-i)nf=i(^-/3.)' 

and for t & E, 



1 



(2.3) ^lim l^(t + .e) =.(*) := _ _ ^^y 

We obtain the following, 

j=0 



(2.4) ^ / W{z)Pl{z)dz -U ^dz = j^rjAn)^. 



where, 



Ja 



I A y{z) 

The Ij can be evaluated by evaluating the residue at infinity, in which case, 
lj = J <g 
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and, 



^g+k 



(^^-i)ro=i(^-c..)(^-/3.) 



1/2 



dz 



which under the substitution z = 1/Z becomes, 

dZ 



ig+k 



I — 



111- 



{l-Z^)UU(^-ajZ){l-(3jZ) 
So in general the integrand of Ig^k has a pole of order fc + 1 at infinity. Hence if, 

-1/2 

fiz) ' - " ' 



then, 
(2.5) 



Ig+k = -2-Ki 



jW(0) 



k\ 



where f^''\Z) denotes the kth derivative of /. Deforming A onto E (2.4) becomes, 



4/i„ = -r]g{n)Ig. 



From (2.5) we have, 

Ig = -27ri, 
and therefore, 

?7g(n) = 2hn- 

To find the remaining rij{n) coefficients we multiply (2.1) by successive powers of z and perform the same 
integration. This leaves us with the following system of equations that we can use to solve for the %(n). 



ro(n) = Igrjgin) 

ri(n) = Ig+ir]g{n) + Igr)g-i{n) 

r2{n) = Ig+2Vg{n) + Ig+iVg-i{n) + IgVgMn) 



Here, 



rgin) = ^Ig+jTJj. 

j=0 



(2.6) rfe(n) = — / z''P^{z)w{z)dz. 
« Je 
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These constants are obtained by iterating the three term recurrence relation k times in (2.6). The recurrence 
relation can also be expressed in matrix form by, 



Lp = zp, 



where, 



L 



( 


1 



















ai 




1 





















b3 


1 



















as 


bi 


1 



















ai 


&5 


1 



























I 


















I 



( Mz) \ 
P2{Z) 



In this case. 



z''p 



L'^p 

OC 



m=0 

oo 



zPn{z) 

Z^Pn{z) = Y.'^L%^P^{z). 
m=0 

From the expression for rfe(n) in (2.6) it is clear that, 
(2.7) rfe(n) = ^[iVi,n- 

To find the ij{n) we use the same method. We know from the proof of theorem (1.3) that, 
^g+i(n) = K-i. 

Multiplying (2.2) by successive powers of z and integrating we have the following system of equations, 

Sl(n) = Ig + 2^g+l{n)+Ig+lig{n)+Ig^g-l{n) 

S2{n) = Ig+3^g+i{n) + Ig+2^g{n) + 7g+i^g_i(n) + Ig^gMn) 



S+1 

Sg{n) = ^Ig+j^j{n). 

j=0 



Here, 



Sk{n) = — ( z^Pn{z)Pn-i{z)w{z)dz. 
« Je 
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In terms of the Jacobi matrix, 



For g = 1 we have, 



Si{x,n) = 2hn{x + bn+i _ + ^) ) n>l 



G2{x,n) = /i„_i ( - ^^^a; + 2a„ - i(a-/3)2 - i ) . n>2 



5i(.x,0) = (x-a) 
G2{x,l) = (x - - a). 

The auxihary polynomials for 5 = 2 can be found in appendix B. 
Non-linear difference equations 

We can use the polynomials studied in this section to derive a pair of non linear difference equations for 
each genus that are satisfied by the recurrence coefficients {a„} and Using the three term recurrence 

relation, 



(2.8) 



= {Z- b„fSn-l{Pz)^n-l{Pz) 

- an-l{z - bn)iSn-2{pz)£n-l{pz) + £n-2ipz)£n-lipz)) + al_i£n~2ipz)£n-2iPz) 

{z - bnYSg{z; n - 1) + a\_2Sg{z] n-2)~ 2a„_i(2 - 6„)Gg+i(2; n - 1) 



Evaluating both sides at z = bn gives, 

Sgibn, n) = af^_-i^Sg{bn, n - 2). 
If we write this out fully for 5 = 1 we have, 

/^^N f, , a + l3\ f, , a + P 

(2.9) an bn+i +bn — = a„_i 6„ + bn-l 



This is our first difference relation for the recurrence coefficients. To get the second we note that, 

m!|(i-a,) = ^(fn(P.)^n-l(Pz)+^n(P.)fn-l(P.)) 

^(((^ - bn)£n-l{Pz) - an-l£n-2{Pz))£n-l{Pz) 

+ {z - bn)£n-l{pz) - an-l£n-2{pz))£n-l{pz)^ 
{z - bn)£n-l{pz)£n-l{pz) ^ n-2 (pz)^n-l (pz) + (pz)fn-l (pz)) 



(2.10) 



{z - bn)Sg{z;n - 1) - a.n-iGg+i{z;n - 1) 

n?=i(^-".) 
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If we evaluate this at z = bn we get , 

Gg+i{bn;n) = -a„_iGg+i(6„; n - 1). 
Writing this out for genus 1, 



(2.11) an + a„-i =2+ 8 + ~ 



This is our second difference equation for the recurrence coefficients. Both equations are valid for n > 3. 

3 Product representation for Pn{x) and Qnix) 

In this section wc use the polynomials derived in section 2 to derive a representation for the generalized 
Chcbyshcv polynomials where the coefficients of the powers of x are given in terms of the branch points and 
the recurrence coefficients. From (1.5) we have 

Pn{x) = ^ (^n(px)+4(Px)) • 

We define the following functions, 

£n{px) 



f+{x;n) :-- 
f-{x;n) := 



Sn-l (Prr) 



Using (1.16) and (1.17) we can write, 

(3 1) f (x-n) = ^"(P'^)^"-i(P'^) = Gg+i{x;n) + hn-iy{x) 
£n-i{px)£n-i{px) Sg{x;n-1) 

Similarly, 

(3 2) f (x-n) = '^"(P^)'^»-i(P^) ^ Gg+i{x;n) - h^^iy^x) 
£n-i{px)£n-M Sg{x;n-1) 

Since, £o{px) = 1, we have, for any x G C 

/ n n \ 

(3.3) Pn{x) = - l^f+{x,J) + l[f.ix,J)\ . 
Similarly, 

(3.4) Q„(^)=M I ]J/^(2;,j)-]J/_(^,j-) 
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(3.3) and (3.4) are valid for any genus g and require only knowledge of the branch points and the recurrence 
coefficients. They do not require any knowledge of the 7 points. If we look at the genus 1 case we find that, 



f±{x\n) = 



^2 _ 



+ 2a„ - |(a -l3f-\± ^(x^ - l){x - a){x - (3) 
2{x + bn+i-^^) 



If we take the limit as a ^ /3 we see that, 



Si{x,n) 2hn{x — a) n> 1 



G2{x, n) hn-ix{x — a), 

and, 

f±{x,n) i(x± - 1) = ie=^*^ a; = cos6' n > 2. 
Note that, 

/±(x,l)^e±'^ 

Thus we can see that (3.3) is a natural generalization of the following formula for the classical Chebyshev 
polynomials, 

{x + V^^T + Vi^^)" 

-l^nKX) — 



2" 

In general, (3.3) and (3.4) are very effective representations of the polynomials for the purpose of extracting 
the polynomial coefficients. The appearance of the theta function of the Riemann surface is still apparent 
in the expression for the recurrence coefficients. One may also wish to study certain special cases of the 
Akhiezer polynomials when the recurrence coefficients become algebraic functions of the branch points a, 
and Pi- In this case (3.1) is particularly useful as one could insert directly the value of the recurrence 
coefficients. In the next section we will encounter some of these special cases. The following formulas follow 
naturally from (3.1): 

and the corresponding formula for the Qn is, 

Letting ai — > (3i in (3.3) we obtain as expected the analogous formula for the Chebyshev case: 

2Tn{x) = xTn-iix) + (a;2 - l)[/„_i(a;). 
Here Tn{x) and Un{x) are monic Chebyshev polynomials of the first and second kind respectively. 
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Now we show how the product representations of (3.1) and (3.2) can be used to re-express the differential 
relations derived in jSj, where Pn{x) and the Qn{x) were shown to satisfy, 

(3.7) P:,{x) = h(x-n)P^{x) + f2{x;n)Qn{x), 
and 

(3.8) Q'„ ix) = hix; n)P„ ix) + h{x; n)Q^{x), 
where, 



and, 

h{x;n) = ip'^{x)f2{x) 

h{x;n) = ^^ + /i(x). 
■ip{x) 

The definitions of all the constants appearing in the above formulas can be found in section 1. From P], 

A(.;n) = ^, 

Pn(X) 

and, 

(x) 

(3.9) f,i:,■n)^^^. 

Using (3.1) we can express fi{x]n) and f2{x;n) in more explicit terms. Recall that, 



Pn{x) 



Sg{x,n) 



Therefore, we can write, 

p'^{x) 1 IS'gix^n) 



(3.10) h{x-n) 



Pn{x) 2 ySg{x,n) ~{x-o^3 
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This is simply a restatement of the earher definition if we recaU that the ji{n) are the zeroes of Sg{x, n). In 
order to re-express f2{x; n), we observe that 



Therefore, 



n 



mn 



From the definitions of f±{x,j) we have, 

'Gg+i{x;j) + h.j^iy{x) 



2i*„(x) (mod 2ti). 



(3.11) 2i^'„(a;) = ^hi 



(mod 2tt). 



Hence, 

n 

(3.12) *<(^)=E 



y'{x)Gg+i{x;j) - y{x)G'gj^^{x] j) 
Sg{x;j)Sg{x;j - 1) 



Substituting the above into (3.12) gives us the differential relations (3.7) and (3.8) completely in terms 
of the recurrence coefficients and the branch points. 

Note that, generally when 7i(fc) € (a;. Pi) k > I, we can evaluate ^'„(ai) by taking the limit x ^ a; in 
(3.11). From (1.13) and (1.14) we see that for x = ai the only non-zero term appearing in the sum in (3.12) 
is the j = 1 term. For j > 1, Gg+i{ai; j) cannot be zero since this would imply, from (1.13) and (1.14) that 
7i(i) = Oil- We have, 

hm 2i^^{x) = lim In ( ^^±lp^{±jM) (mod 2^). 
x^a, x^a, \Gg+i[x;l) - y[x) J 

We know that, 

g 

Gg+i{x; 1) = (a; - h) J|(a; - aj) 
and therefore, 

7r 

lim 5',i(a;) = — (mod n). 

x^ai 2 



3b. Discriminants 

In this section we derive an expression for the discriminant when g = I. Once we have knowledge of 
the differential relations satisfied by a sequence of polynomials one can in general use this information to 
say something about the discriminant. The discriminant is useful when we consider certain electrostatic 
problems regarding the zeroes of the polynomials. See and ^H] for more details. The expressions we 
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derive in this section are restricted to the cases when the polynomials do not have zeroes at the branch 
points. If we write, 

n 

j=l 

then the discriminant is defined to be the following, 
(3.13) D[P^{x)]:^ \[ (x,,„ - a;,- „f . 

\<i<j<n 

Since 

Pn{x) 



it is easily verified that, 

n 

{-l)-^\{Pi{x,,n)=D[Pn{x)]. 
J=l 

From (3.7), assuming that Pn{x) has no zeroes at the branch points, we have, 

PniXjjn) ~ fliXj^n^QniXj^n)- 

To evaluate Qn{xj,n), recall that, 

Pn{x) = Pn{x) COS ^n{x), 

and 

Qn{x) = i^|J{x)pn{x) sin <i/ri{x). 
Since we have, 

1 Qn{x) 
Wx)TJx)=''''^-^'^^ 

it must be that, 

tan*„(a;j,„) = ±oo, 
and therefore if Xj^n is a zero of Pn{x) then we do have, 

cos*„(xj,„) = 0. 

Hence, 

Sm'^ n{xj,n) = ±1- 
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Therefore, it is always true that, 



\Qn(Xj.n)\ = \ilp{Xj,n)Pn{Xj^n)\- 

The sign depends on the value of n and j. To show this dependence, we note that the following quantity. 



2K 



Xj,n - li{n) 



is positive. Since theorem (7.1) in [S| states that any zero of Pn{x) on the interval (a/,/3/) must lie in (a;, 7;) 
and since the zeroes of P„(x) and Qn(x) interlace, we must have 



and, 



sin^'2„(a;j,2n) = (-1)^ 



3in^'2„-i(xj-2„-i) = (-1^+^ 



Therefore, 



and, 



D[P^{x)] = {-l)^X{PL{x,,n) 

n 



"V (4"-i)nf=i(^.,n-A) 



f2{x3,n) 



(3.14) 



= (-l)t(2/i„)^ 



n?=l^n(7,(n)) 



^n(im.(-i)m=iJ^n(/3, 



W_h{Xj,n)- 



Now in this most general form we can not yet evaluate the discriminant for general g since we do not 
know how to evaluate the product of the /2(a;j,„). However, examining the form of f2{x] n) we see that it is 
a rational function. In certain special cases we may know enough about this function to factor it in which 
case we can express the unknown quantity in (3.17) as Pn{x) evaluated at the zeroes and poles of f2{x;n). 
For g = 1 we can always do this, in which case. 



D[Pr,{x)] = (-l)t(2/i„)t 



P„(l)P„(-l)P„(/3) f \ 



Y[f2{Xj,n)- 



Now, 



f2(x;n) 



n{x — r_|_(n))(x — r_(n)) 
{x — j{n)){x — a) 
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where, 



r± (n) 



(co/n + 7) ± V(co/»^-7)^ +2yiA 



Therefore, we have for the square of the discriminant 

(-l)"P,?(r+)i^2(^_) 



(3.15) D'[Pn{x)] = (2n2/i„)" 



p2(«)^n(l)^n(-l)P«(/3)P„(7(n)) 



The product Pn{P)Pn{ln) is always positive since we know from 9 that Pn{x) cannot have a zero in the 
interval (7(71), A similar calculation for the polynomials of the second kind gives. 



(3.16) D''[Q^{x)] = (2n2/,,„)^ 



Ql{T+{n))Ql{r.{n)) 



Ql{l)Ql{~l)Ql{P)Qn{l{n))Q.r,{a) 



Again, in this case we require that Qn{x) have no zeroes at the branch points. 

These discriminants are valid for genus 1. Indeed we know explicitly the values of 7(n) and r±{n) for 
g = 1. It is clear from the derivations in this section, that the equivalent expressions for general genus g are, 

(3.17) D^mx)] = {2n^KT ("l)" -li ^'n^C'^, W) 



Pn{l)Pn{-l) Y{U^PnmPn{n3{n))P^A^,)] ' 



and, 



2, \n^l 



(3.18) D'[Q^{x)] = {2n'hn) 



Q?.(l)Q^(-l)nf=i[Q^(/3.)Qn(7,("))Qn(a,)]^ 



where the {rj{n)\ are the 2g roots of the polynomial. 



W_{x--fj{n)) 



9-1 9 
nx'' — > Cj{n)x^ > 



2 ^^x--ij{n) 



4 Connection with Polynomial mappings 

The theta function expressions for the recurrence coefficients a„ and hn of the generalized Chebyshev polyno- 
mials are given in (1.5) and (1.6). Consider for the moment two different values of the recurrence coefficients, 
for example, 6„ and fon+^Si where n and ^ are both integers. If ^ is chosen such that .SB G then from the 
periodicity properties of the theta functions we have, 



5ri+i? — On- 



From the expression for the a„ coefficients, we will also have. 



an+^ = a„ n> 2. 
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When the entries of the vector B are all rational numbers, we can always find an integer K such that G . 
If we take the smallest integer that satisfies this requirement we will then have the smallest integer over which 
the recurrence coefficients themselves repeat, that is the recurrence coefficients are periodic with period ^. 
If at least one of the entries is an irrational number we are unable to find an integer R that satisfies the 
periodicity requirement. 

We will study the effect that this periodicity has on our previous constructions and use this to explain 
the connection between the generalized Chebyshev polynomials and non-trivial polynomial mappings |10| . 
Throughout this section is a given fixed integer. First consider Pnii{x). From (1.5), and with the periodicity 
condition, we see that the 9 factors cancel in equation (1.5), therefore 

(4.1) P^^{x) = ^(e"^("(P^)) + e-«^(^2(P.)))^ 
hence, 

(4.2) P^(x) = ^(e-s(f^(p^)) + e-«(f^(P^))). 
Here, 

(4.3) A„^ = 2e-"^'^''. 
Therefore, 



(4.4) = A'^t/^^^''^ 



V A, 

where Tn{x) is the monic Chebyshev polynomial. From the expression for hn in (1.10) we see that. 

It follows that, 
A^ = 
and therefore, 

(4.5) P„^(x) = (2/1^) tT„ ' ^-'^^'^^ 



The corresponding formula for the polynomials of the second kind is, 



(4.6) Q„^(a;) = (2/1^) — Q^(x)C/„ , 
where Unix) is the Chebyshev polynomial of the second kind (degree n — 1). Recall that, 

^n(px)^n(pi:) = TTn T- 

nLi(a;-ai) 
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When ^ e Zs, we have, 



^„(px)4(p.)=4e-2»>^''. 

Therefore, since the 7j(n) are the zeroes of Sg{x; n) we must have 7j(^) = aj. It wiU be shown later in this 
section that the maximum value of g permitted for the polynomials is ^ — 1. Therefore, if we put g = ^—1 
and n = ^ in (1.13) we see immediately that, 

Qa{x) = JJ(a;-aj). 



Hence, 



Psi{x) 



In the language of (3.3), (4.5) becomes, 

p , , (n^^i/+(x,j))" + (njLi/-(^,j))" 

and (4.6) becomes, 

^ , , , (nf=i/+(^,j)r-(njli/-(^,i)r 

Using the theta function expressions, we can verify the following, 

f+{x,R+l) = = ^^/+(x,l). 

Similarly we can deduce that, 

Afi+i ^ 1 
A^Ai 2' 

and therefore, 

f±{x,^+l) = y±{x,l). 
An analogous calculation shows that, 

f±{x,^ + j)=f±{x,j), l<j<R, 

as we expect from looking at the explicit structure of the f±{x,j). 

If we know the period ^ then we can use the Pg(a;) with (4.5) to generate every polynomial in the 
sequence. To obtain the intermediate polynomials PnA+j for j = l.M— 1 we use the fact that 
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and therefore, 



1 / 1 

(4.7) Pn^+jix) = - yPj{x)PnR{x) + ^i^Qj{x)QnM{x] 



l<j<A, n > 1 



Similarly, 



(4.8) QnM+j{x) = -{Pj{x)QnA{x) + Qj{x)PnA{x)), 1 < j < ^, n > 1 

So if we know the first K polynomials of the first and second kind we can generate all of the remaining 

polynomials using (4.5), (4.6), (4.7) and (4.8). 

From the theory of the classical Chebyshev polynomials we have, 
T„(a;)T„(a;) = T„+„(a;) + ^T„_„(a;), m > n ^ 



Similarly, 



Um{x)Tn{x) 



Um+n{x) + ^Um-n{x), m>n^O 
Um+n{x) - oL:Un-mix), n > m ^ 



Analogous formulas to the ones above do not hold for all of the generalized Chebyshev polynomials. However, 
PnA play a similar role to the classical Chebyshev polynomials. 



Lemma 1 For any m,n G Z+ we have 

(4.9) Pn{x)PmM{x) = 



Pn+m^ix) + i^Pn-mSi{x) U > mA ^ 



Pn+mM{x) + ^PmSi-n{x) toJ? > 71 ^ 



Proof : 

From section 1, 

Pn{x) = e~ 



■nxo _ 



t?(u+) 



i?(u+ - nB) 



^nn(p.) ^K-nB) ^ ^_„o(p^) ^^(ux + nB) 



1?(Ux) 



i?(u.) 



and, 



Pm^ix) = e~"^^^° gm.Sn(p»,) _|_ g-m.sn(p»,) 



Therefore, 



Pn{x)PmA{x) 



??(u+) 



i?(u+ - nB) 



'^(n+mj^)n(p.) -^("^ - "B) ^ ^_(„+^;^)o(p^) t^(u^ +nB) 



1?(Uc. 



^(Ux) 



1?(Uc.) 



t5(U:, 



From the periodicity of the thcta functions wc may add the vector ±m.SB to any of the arguments of the 
theta functions without changing the equality. We make this change in the argument of the theta functions 
in such a way that we are left with (4.9). □ 
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Similarly, for the polynomials of the second kind, we have 

Qn{x)PmA{x) = 



Qn+mSi{x) + ^^Qn-mfiix) U > TTlR ^ 

Qn+mA{x) - ^QmA-n{x) m^>n^O 



Combining (4.9) and (4.7) we obtain, 

P„^_,(a=) = 1 (p„.(a=)P,(a=) - ^^^^^^g^') j = l,-,^- 1. 



Similarly, 



QnA-j{x) = ■^{QnA{x)Pj{x) - Pn^{x)Qj{x)). 



From (1.8) it follows from the periodicity properties of the theta functions that. 
The bn also have a similar behaviour, and one can show that 

However, this is not obvious from the form of the recurrence coefficients given in (1.9). To prove the above 
equality we start with the three term recurrence relation (1.4) with n = A — j + 1 and a; = 1, so that 

Ps_,+i(l) = Pft-,+2(l) + 6^-,+2Pg-,+i(l) + a^_,+iPs_,(l). 

From the expression for PnA-j {x) we have. 



hj 

Substituting this into the three term recurrence relation we obtain, 

P,_i(l) = Pj{l) + 6^_,+2P,-i(l) + a,_iP,_2(l). 
Comparing this with (1.4) at a; = 1 and n = j — 1 we see that 6j?-j+2 — bj. 

Theorem 4.1 Given a set E such that the recurrence coefficients have period ^, then the polynomial 
Ps(a;)/Aj^ maps each sub-interval of E onto [—1,1]. 

Proof: Take an arbitrary polynomial M^{x) of degree ^ with leading coefficient l^ where all of the zeroes 
lie in (—1, 1). Furthermore, 

\M^{xi)\ > 1 



26 



where the the stationary points, and, 



M^(l) = 1 
M^(-l) = (-1)^. 

Then the inverse image of [—1,1] under Mf^ is a set J of at most disjoint intervals of the form J := 
[—1,0;] ufZi [l3j,aj+i] U [/3j?_i,l], with —1 < aj < (3j < aj+i < 1. Let arix) be the weight function 
associated with the Chebyshev polynomials, 

aT[x) = -- 



TT \/l - 

Following ^U] we define a function W{x) by, 

aA{x) = W{x)aT{M_^ix)). 



(4.10) 



nf=i^(a;-aj) 



(i-^')nf=7(^-/3.) 



Now, M^{x) has the property that it maps all of the branch points into —1 or 1, and therefore, 

A-l 

Mlix) - 1 = iiix' - 1) n - - 

i=i 

Since Msi{x) is a polynomial, the above equation imposes constraints on the branch points. Substituting for 
M^{x) in (4.10) we have, 

Si-l 

W{x)^lsiY[ix-a,). 
i=i 

Now 



(4.11) 



W{z) 



Wi{x) 



Mj,{z)-x Utii^-M-\x)) f-[z-M-\x) 
where the M~^{x) are a complete set of inverse branches of Msi{x). Now, 



W{z) = (M^(z)-a;)^ 



Wi{x) 



and 



-{z-M^ \x) 



W{M-\x)) = w,{x) Y.{M-\x) - Mj\x)l 



and therefore. 



W{M^\x)) 
M'{M-\x)y 
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Also, by comparing coefHcients of ^ in (4.11), 
J? 

^Wi{x) = 1. 

From the orthogonality relationship of the Chebyshev polynomials we have, 

hl^^Snm = J Tn{z)Tra{z)(JT{z)dz 

= jT^{Ms^{x))T^{Msi{x))(JA{x)dx, 
where the last line follows from (4.10). Using the general method in JUj theorem 1, it can be shown that 

x-'Tn{Msi{x))aA(x)dx — j < nM.. 
Therefore, 

T„(M^(x)) = CnPndx) n>0, 
where c„ is some constant. So for n = 1, 

Msi{x)^ciPji{x). 
Now, from (4.5), 

and therefore ci = and we have proven that ^-^^ maps J into [—1,1] A -fold. Also, since for the 
polynomials orthogonal on the set J we have the shown (4.5) to hold, we can conclude that the recurrence 
coefficients are periodic with period ^ and therefore J is in fact the set E referred to in the statement of the 
theorem. □ 

Remark 1. This theorem has appeared in a more general setting, see |14| and |15j where necessary and 
sufficient conditions were given for the existence of an extremal polynomial on a certain set E, which corre- 
sponds to our mapping polynomial. In general the orthogonal polynomials on these sets have asymptotically 
periodic recurrence coefficients. These are found to be orthogonal with respect to a more general class of 
weight functions similar to w{x) of the form. 
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where Trji-i{x) is a polynomial of degree ^—1 chosen such that the weight function is positive on the interior 
of the set E. In equation (1.12) of 15 the extremal polynomial of degree n, Tn{z), is where n, the 

asymptotic period is identified with our 

As a consequence of this theorem, when the recurrence coefficients have period K the set E has at most 
R sub intervals. Recall that g is the number of gaps in [—1,1]. Therefore, the maximum value of g for this 
set is — 1. 

In figure 1 we show an example of a mapping polynomial for A = 3. 




Figure 1: General Form of Af3(a;) compared with the Chebyshev polynomial of degree 3. The branch points 
are where the polynomial crosses the lines y — —1 and y — I. 

We now consider some examples of sequences of orthogonal polynomials generated via the polynomial 
mappings P^{x). We choose small values of R and investigate. M. = 1 corresponds to the Chebyshev 
polynomials as the recurrence coefficients are constant. 

For R — 2, the mapping polynomial, 

M2{x) = -^P2{x) = -^{x^ +ax + 62), 
satisfies the following conditions 

This implies that M2{x) has the following form 
with 



A2 = — — = — 
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Therefore, 



1 — \ 1 — / 

with the condition on the branch points, a = — /3. The condition above, tells us which sets E arise from a 
polynomial mapping of degree 2. These are the same sets that give rise to periodic recurrence coefficients 
with period 2. To obtain the recurrence coefficients for these polynomials, we can substitute f3 — —a into 
the general genus 1 recurrence coefficients. The general expressions for the recurrence coefficients can be 
obtained in one of two ways. Appendix A shows how we can manipulate the theta function expressions for 
the recurrence coefficients into expressions involving the Jacobian elliptic functions. We are then able to 
invert the functions and in principle obtain the recurrence coefficients for any n in terms of the branch points 
only. In this case we only need to find the first two since the others repeat. The second way to find the 
recurrence coefficients is to use the non- linear difference equations derived in section 2. Once we have the 
initial conditions we can iterate the non-linear difference equations to find the remaining coefficients. This 
method is more practical when we are considering higher periodicity. In any case, the recurrence coefficients 
are found to be, 

(4.12) a„ = n>2. 
With these coefficients, our polynomials become, 

and 

Also the auxiliary polynomials become, 

Si{x,n) = 2/t„(a; + (-1 
G2{x, n) — hn^i{x + a){x — a). 

From (3.10) we have, 

fiix;n) = i ( — , ^ —] ■ 

2 \x + (— l)"+^a x — a J 

To evaluate f2{x] n) we must know the co{n) and ^{n). For genus 1 polynomials we have from [H], 

a + f3 j{n)-a 
Coin) = pi{n) + + ^ • 

where pi{n) is the coefficient of in Pn{x). In this case, 

pi(2n) = 
pi(2n — 1) = a. 
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Therefore, 



co(n) = 0. 
Since 7(n) is the zero of Sg{x, n), 
7(n) = (-l)"a, 

and, 



. , . nx 
f2{x;n) = 



X — a 



For even n, the discriminant becomes (n = 2m), 

D[P2mix)] = (-l)'"2-(2'"-l)(2™-2)-m(2TO)2™(l - a2)™(2'"-l)T„ ( 

— 1 

and for odd n, {n = 2m + 1) 

D[P2m+l{x)] = (-l)'"2-*'"'+'"-l(2m + 1)2™(1 - 02)^(2™+!) 

a2 + l\ 2 /a2_^i 



a2 — 1 y 1 — 0:2 — 1 

The discriminant for the case when n is odd cannot be obtained directly from (3.15) since the P2n+i{x) have 
zeroes at x = —a. The problem is that of evaluating P2n+i{~'^)- general evaluating the derivative at one 
of the branch points is not straightforward. However, in this case it can be verified that, 

PL+ii-ct) = 2/2(-a;2n + l)Q2n+i(-a). 
whereas, for the zeroes that do not lie on the branch points we have, 

P2n+l{^j,'2n+l) = /2(a;j,2n+l ; 2n + l)Q2n+l{Xj,2n+l) 

and the factor of 2 docs not appear. It is for this reason that the expression given in (3.15) remains valid in 
this case if wc insert a factor of 2. 

The polynomials in this example can also be shown to satisfy the following differential equation, 
x{x-^ — a-^){l — x-^) [1 — x-^){x-^ — a'^) 



and. 



„ a;^ - 2aa;(l - a;2) + n'^x'^ - 2xa + {v? + l)x^a - o? „ _ „ 

" ~ a;(a;2-a2)(l-a;2) x{l - x-^){x - a){x + " ~ 

for odd n. 



For the next example we choose .ft = 3. Therefore our polynomials can be described in terms of a 
polynomial mapping of degree 3. The highest genus permitted in this case is 2. This corresponds to three 



31 



disjoint intervals. Genus one cases arise if we close one of the gaps. In the last example we saw that the 
only sets that gave rise to a period of 2 were symmetric around the origin, that is a = — /?. The aim of 
this example is to determine which sets give rise to recurrence coefficients with period 3. We determine the 
recurrence coefficients, the 7 points and under a certain condition, the Cj{n). 

We have a mapping polynomial M3{x) with leading coefficient ^ that satisfies the following equations, 
M3(±l) = ±1 

Mi{a2) = M3(/32) = -1. 
Solving these six equations we have, 

^3(^) = £+«^'+(l-^)x-a, 



with. 



(4.13) a = 



+ l)(ai + 1) (/32-l)(a2-l)' 
and, 

(4.14) A3 = + l)(ai + 1) = {P2 - l){a2 - 1). 

Since we can characterize the mapping polynomial in terms of two branch points, we have a system with 
two unknowns. We describe the system in terms of ai and f3i. The form of the mapping polynomial restricts 
ai and (3i to certain intervals on [—1, 1]. If we eliminate a2 from (4.13) and (4.14) we get, 



^ ai + /3i + 2 ± V(/3i - ai? - 4(ai + /?i + 1) 
02 = ^ 

hence, 

(4.15) (ai-/3i)2-4(ai+/3i + l)>0. 
Therefore we have. 



ai < /3i + 2 - 4 



2 

This inequality, coupled with the inequality, 

implies that ai lies on the interval [—1, —5]- Therefore /3i lies on the interval [ai, /3i^max] where. 



01,max = 2 + ai - 4 



ai + 1 
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The value of Pi^max follows from (4.15) with the equality sign. With these conditions on ai and /3i we have, 

(4.16) a2 = 1 + - l^/iPi - ai)2 - 4(ai + /?i + 1) 
and, 

(4.17) P2 = l + + l^/iPi - o^iY - 4(ai + A + !)• 

For g = 2 we evaluate the first three recurrence coefficients in the most general case in terms of ai, a2, Pi 
and (^2, using (1.13) and (1.14). Then we substitute for a2 and (32 using (4.16) and (4.17) and find. 



i(l + ai) (2 - ai + (h + v/(/?i-«i)'-4(l + ai+/3i)) 
^ (ai - /3i - 2 + v'(/3i-ai)'-4(l + ai+/3i))' 
(4.18) as = J(l + ai) (2 - ai + /Ji + - ai? - 4(1 + ai + /3i)) , 



Oi 
02 



and. 



61 = ^ - ai + V(/3i - ai)2 - 4(1 + ai + /3i) 

62 = ^ (2 + 3ai + /?i - v'(/3i - ai)^ - 4(1 + ai + /3i)) 
(4.19) 63 = ^ (2 + 3ai + /3i - - "i)^ - 4(1 + ai + /3i)) 



We also find that, for n > 2, 



7i,2("') = ^il + ai+Pi-bn+i 



T \Jl- 2(ai + /3i) + (ai - /3i)2 - 4(a„ + a„+i) + 26„+i(l + ai + /3i - ^^n+i) j • 

Expressions for the polynomials are easily obtained from our equations in the previous sections. Expressions 
for the discriminant and the differential equation can in principle be obtained using the constants that we 
have derived in this section. 

A highly symmetric case corresponding to g = 2 arises when (3i + ai = —1. With this, the a defined by 
(4.13) is 0, which implies that 

a2+ P2 = 1- 
Substituting these into (4.14), gives, 

P2 = l+Pl. 
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This case corresponds to the first gap being symmetric about the point —-^ and the second gap being 
symmetric about the point 5. In this case the square roots appearing in the recurrence coefficients vanish 
leaving (with a = a{), 

a\ = — 2a(l + a) 
1 

a2 = - 

03 = + a) 

and 

61 = -(l + 2a) 
h = ^ + a 

Also, 

G3{x,n) = hn-i{x^ + (2a„ -1- a- a'^)x + 2an{bn + bn+i)) 

and 

S2{x, n) = 2hn{x'^ + bn+ix + a„ + a„+i + - I- a - a^). 

This gives 

-hn+i T J bl^i - 4(a„ + a„+i + 6^+1 - 1 - a - a^) 
7i,2(ri) = . 

The first three polynomials are, 

Pi (a;) = x + 2a+l 

P^{x) = x2 + (a+i)x-i 
P^{x) = x^ - {1 + a + a^)x. 

From (3.9), (3.10) and (3.12) we can calculate fi{x;n) and f2{x;n). 

Returning to the previous non-symmetric genus 2 example, we add the restriction that one of the gaps 
is closed. We take the case where a2 = ^2- In this situation (4.13) and (4.14) simplify to 

^ -(/3i+Qi + l) ^ -(2a2-l) 
(/3i + l)(ai + 1) (a2-l)2 ' 

and 

A3 = (/?i + l)(ai + l) = (a2-l)^. 
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Therefore, eliminating a-i from the above we obtain, 

(/3i + l)(ai + 1) = 
with the solutions, 

/3i = ai + 2 ± 4^ 
Since a\ < (3i we must have 

(4.20) /3i = ai + 2 - 4 



a.i + 1 



ai + 1 



This relation governs the variation of /3i as ai varies over the interval [—1, —5]. It describes a contour in 
the (ai,/3i) plane. Figure (2) shows a graphical representation of this contour with similar curves for higher 
values of ^. Had we have chosen the condition ai = (3i for the other g = 1 case we would have obtained, 



(4.21) /?2 = -2 + Q2+4 



1 — 0:2 



where 02 varies over [—1,^]. Using the method outlined in the appendix we are able to compute the 
recurrence coefficients for 5 = 1 explicitly. From now on we refer to ai and (3i as just a and (5 respectively. 
Upon substituting for (3 given by (4.21) into the general recurrence coefficients we obtain the following, 



ai = 2(a + l) 1 



l + a 



a + 3 ^ /a + 1 
«2 = ^-2) 



as = (a + 1) 1 



l + a 



and. 



61 = 1-2 

62 = (a + 1) 

63 = (a + 1) 



l + a 



l + a 



l + a 



The recurrence coefficients repeat in blocks of three however we must remember that azn+i = n>\. 
Furthermore, 



co(3n) = n I a — 2 



l + a 



co(3n+l) = M(l + 2n)a- (l + 4n)i 



co(3n + 2) = M(l + 2n)a- (3 + 4n)/ 
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and, 



7(3n) = a. 
7(3n+l) = 

7(3n + 2) = 

These constants, listed above, allow one to calculate the differential equations and the discriminant given in 
the previous section. 

General A 

For notational convenience we now label as the sets E that give rise to period ^ recurrence coefficients. 
Equivalently these are the sets for which the minimum degree of polynomial mapping from onto [—1, 1] 
is ^. Prom our earlier examples, 

£1 = [-1,1] 

E2 = [-l,a]U [-a,l] ae[-l,0] 

E^ = [-l,a]U[/3,u{a,/3)]U[v{a,P),l] 

where, 

u{a,p) = l + ^-^V(/3-")'-4(a + /3+l) 

v{a,f3) = l + ^±l + ^^{i3-ar-A{a + (3+l). 

Each set i?^ is parametrized in terms of .ft — 1 variables. As R increases it becomes intractable to solve 
explicitly for the relationships between the branch points. However, the branch points are constrained on 
certain curves in [—1, l]'^^ where the co-ordinate axes of the 2g cube label the position of the branch points. 

In what follows, we describe the constraints on the mapping polynomial. This will provide a set of 
equations that define the curves of constant periodicity in the 2g cube. We also show how to determine the 
smallest R for a given configuration of the branch points. 

For general the conditions imposed on the mapping polynomial, 
M^{x) = a^x^ + aR-ix^~^ + ... + a-^x + ao, 

are 

M^(-l) = (-1)^ 
M^(l) = 1 

(4.22) M^iaj) = M^{f3j) = (-1)^+^' j = 1, ...,^-1. 

Therefore .ft + 1 of these will determine the coefficients in terms of .ft + 1 branch points. Using the first two 
equations of (4.22) and any .ft — 1 of the remaining ones we can solve for the coefficients of the mapping 
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polynomial in terms of 1, —1 and ^ — 1 free parameters. The other R — 1 branch points are related to the 
free parameters by the remaining equations of (4.22). These equations determine a surface in the 2g cube. 
The surface obtained in this way will correspond to maximum g = M. — 1. Closing one of the gaps, that is 
setting aj = /3j for some j < R — 1, will determine surfaces corresponding to that value of M. for smaller 
values of g. 

If we close all of the gaps by setting aj = 13 j for all j = 1....^ — 1, then the mapping polynomial becomes, 
M^{x) = 2^-^T^{x). 
If we denote by Xj the point where aj = f3j, then we must have 
2^-'T^{xj) = ±1 

from the properties of the mapping polynomial. Since Tn{x) = 2^^" cos{n0), x = cos{9) we see that the set 
of points {xj}fzl are the ^—1 stationary points of the Chebyshev polynomial of degree ^. Therefore, each 
gap i.e. {ai,l3i)fz^ can be said to be 'centred' (not necessarily symmetrically) around one of the stationary 
points of the Chebyshev polynomial of degree ^. The most simple example of this is the ^=2 case where 
the gap is centred around the point x = 0. In the ^ = i, g = 2 case, the two gaps are centred around x = ^ 
and X = The stationary points of T4{x) are x = a; = 0, and x ~ So for a polynomial mapping 

of degree 4 we have three distinct 5 = 1 cases corresponding to one gap which is opened around one of the 
stationary points. 

Consider now the ^=2 mapping polynomial, which maps [—1, a] U [—a, 1] to [—1, 1]. Prom (4.4), 

Since, A2 = 2A4, and from the fact that the classical Chebyshev polynomials map [—1, 1] to [—1, 1] we see 
that, -^^^ also maps [— 1,Q;] U [— a, 1] to [—1,1]. Therefore, it is this configuration of the branch points 
which corresponds to the R = A, g = 1 case where the gap is centered around a; = 0. Indeed for any mapping 
polynomial of even degree, one of the g = 1 cases corresponds to this configuration, since a Chebyshev 
polynomial of even degree has a stationary point at zero. 

The same reasoning as above applies for general K that is, if is a mapping polynomial for a 

certain configuration of the branch points, then so are the polynomials Therefore amongst all of 

the configurations that correspond to a mapping of degree nR, there can be found those that correspond to 

degree jR, where j divides n. 

If xi is a stationary point of T^{x) := 2^~^T^{x), then 
and 

which indicates that xi is also a stationary point of Tjj^(x). For example Ti2{x) has amongst its stationary 
points, the stationary points of Tq{x), T4{x), Ts{x) and T2{x). 
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So as an example, a mapping polynomial of degree 12 will give rise to eleven g = 1, single parameter 
cases of the generalized Chebyshev polynomials. Six of these will already be the same as those that arise 
from a mapping of degree 6 and four will be the same as those that arise from a mapping of degree 4. All of 
these cases cover those that arise from mappings of degree 2 and 3. Therefore a mapping of degree 12 will 
give rise to two 5=1 cases that cannot be generated from polynomial mappings of lower degree. 

The upshot of the above discussion is that we can uniquely label each surface in the 2g cube by a set of 
points corresponding to the stationary points of the Chebyshev polynomials of any degree. We do this in 
the following way. 

If we denote by s; the ordered set of all stationary points of the Chebyshev polynomial of degree 1, then 
we can label a surface by taking any subset of the s;. The number of elements in a subset is the number of 
gaps g in the interval [—1, 1]. For example, 



S2 = {0} 



«5 = |-i(v^+l),-i(y5-l),i(V5-l),i(y5 + l)|. 



Therefore, the set, 

1 1 

71' 71 

corresponds to the R = 4 g = 2 surface where both gaps are centred around the two points — ^ and 
This is a two parameter curve. Those curves shown in figure 2 can be labelled by one number which is a 
stationary point of a Chebyshev polynomial. It is clear from the discussion so far that the number J\f given 

by, 
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is the number of surfaces in the 2g cube associated with the integer 

So far, we have discussed in general terms how to describe the surfaces corresponding to a particular 
value of R. On the other hand, the problem of finding R from a given set E appears to be a very difficult 
numerical one. The condition for periodicity of the recurrence coefficients is 

(4.23) e \ {QS}. 



where, 

1 



^ 2m ' 
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If — Bj are rational then R is the least common multiple of the denominators of these national numbers. For 
example, if 




then ^ = 30. However, if at least one of the components of the vector B is an irrational number then it 
is impossible to find an integer K that satisfies condition (4.25). These correspond to the cases where the 
recurrence coefficients never repeat themselves or equivalently, the sets for which there is no polynomial 
mapping from this set onto [—1, 1]. 

Note that in general the values of Bj are subject to certain conditions. In 0|, section 10, a physical 
interpretation was given; the values — Bj were seen to be the proportion of the overall charge lying on 

the interval Ofj], where Pq := —1, in the equilibrium charge density for the set E, when the total 

charge was normalized to 1, that is, 

a{x)dx = 1. 

Since the total charge on all of the intervals must add up to 1 the following inequality is satisfied by the 
components of B, 

0<^-B,<l. 



To finish this section we refer to figure 2. Recall that each point in the plane represents a configuration 
of the branch points for g — I and therefore a sequence of orthogonal polynomials. Using our notation, if 
{xi} represents a curve in this plane with integer ^ then {—Xi} represents a different curve but with the 
same integer ^. We ask the question "Is there a relationship between the polynomials associated with the 
curves {xi} and {— Xi}?" A similar question can also be asked for g > 1. The relation between these cases is 
given in the most general form by theorem 1.2. However, since the recurrence coefficients are periodic, the 
elements appearing in the determinant simplify: 

P„^+, (A) = ip„K(A)^, (A) 1 < J < J? - 1 

and 

where, 

Pniiiai) = PM) = ^^^„!, . 

As an example we show the relation between the polynomials associated with the curves {—5} and {^} which 
are the two R = i g — \ curves. Using the notation of theorem 1.2, we denote the polynomials associated 
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Figure 2: Contours of constant periodicity in the {a, (3) plane for genus 1. Here a = (-1 + \/5)/4, b = 1/2, 
c = l/\/2 and d = {1 + \/5)/4. This is not a precise plot, it is representative of the relative position of the 
contours. 

with the set { — 5} as Pn{x), the polynomials associated with the set {|} as Pn{x), and the respective 
mapping polynomials as Ms{x) and M3{x). It is easy to verify that, 

Msix) = {-l)M3{-x). 

Consequently, 

Psnix) = i-lf^Psni-x). 

However, for the intermediate polynomials we must use the expressions given in theorem 1.2, which give. 



hn+lix) = 



(-1) 



3n 



a + 2-4 



'a±T\ 

2 ; 



P3n+2{-x) - 
+ ( (a + 3) 



a + 1 



-l]P3n+l{-x) 



a+1 



2{a + l)\Psn{-x) 
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and, 



/I ^3r^+l 



4b. General form of the polynomials 

In this section we show some plots of the polynomials for small ^ and prove a theorem regarding the 
zeroes. The Chebyshev polynomials T„(x) is bounded between 1 and —1. These are the mapping polynomials 
from [—1, 1] to [—1, 1]. As a consequence of (4.5) we see that |Pnj?(a;)| < 1 for x £ E, where, 

Outside of the set we no longer have this bound. Therefore we would expect the polynomials to appear 
graphically as in figure 3. In figures 4, 5 and 6 we plot some of the intermediate polynomials PnSi+j{x) for 
various values of j and .R, where, 

In E they are bounded above and below by some function. It is clear from the representation (1.18) that 
the form of the envelope is ±p„j5+j(x), for x € E where. 



PnSi+j{x) — 



I Sg{x, + j) 



(A.s)"y2H7Y ULii^~ai)' 

V n?=i(^-".) ■ 

Note that the quantity under the square root is strictly positive for x € E. 

We already know something about the zeros of the generalized chebyshev polynomials in the most general 
case from In the polynomial mapping case we can say a little more about the zeros. 

Theorem 4.2 All of the zeros of PnSi{x) lie inside E, and in each of the R sub-intervals that make up E lie 
n zeros. 

Proof: All of the zeros of the Chebyshev polynomials are contained in (—1,1). From (4.5), and since 
^-^^ G [^li 1] if only if x £ E, we see that all of the zeros of Pnnix) lie inside E and each subinterval 
Ei contains n zeros. 

The zeros of Qnnix) behave similarly except in this case there are n — 1 zeros in each interval Ei and 
there are also zeros at every aj, j — 1, 1. This is easily verified by checking the form of (4.6). From 

(4.6) the n — 1 zeros in each Ei come for the zeros of Un ('^'J')- each interval Ei, these zeros interlace 
with the zeros of P„^(a;). 
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Figure 3: R=2, j=0 n=1..8, a = -0.5 



Theorem 4.3 Let {xi^n} denote the zeros of Pn{x) and {yi,n} denote the zeros of Qn{x). Consider the 
following two sets, 

A := {{a;i,„^}, {xjj}} 

B := {-l,l,/3i..../3i?-i,{2/ij},{2/i,„j?}}\{ai...a^_i}. 

The following three statements are true : 

Between, any zero of Pj{x), and any element of B which is larger than this zero, there lies at least one 
zero ofPnA+j{x)- 

Between any zero of Qj{x), and any element of A which is larger than this zero there lies at least one 
zero ofPnA+j{x). 

If Pj {x) has a zero at any of the points in B then PnA+j {x) has a zero at the same point. 
Proof: From (4.7), Xi^nH+j solves the following equation, 

PnM{x)Pj{x) = ~-i^j-^Qr,Si{x)Qj{x). 

Using (4.5) and (4.6), we see that this is equivalent to, 

(4.24) A^Tr. ^^•(^) = (^) - n - 

For clarity, we illustrate a proof of the case of = 3 and n = 2. The zeroes of C{x) and D{x) where, 

C(.):=A,r„(i|i)) 
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-4 - 



Figure 5: The symmetric case for ^=3, j=0, n=17, a = —0.7 



and, 

are indicated in figure 7. The crosses denote the zeros of C{x) and the dots denote the zeros of D{x). We 
see that they interlace. In order to plot the remaining zeros of (4.24) we must add j crosses and j — 1 dots, 
corresponding to the zeros of Pj{x) and Qj{x) respectively. We can distribute these without any restriction 
other than they must also interlace with each other. The result is a sequence of crosses and dots that 
alternate except for j occurrences of two adjacent crosses and j — 1 occurrences of two adjacent dots. Where 
we see two adjacent crosses, one of them is a zero of Pj{x). Similarly with two dots, one of them is a zero of 
Qj{x). As an example we plot where the zeroes of both sides of (4.24) might lie, for .ft = 3, n = 2 and j = 2 
in figure 8. 

In order to graph the polynomials that pass through the zeros we need an initial condition. The condition 
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Figure 7: The zeroes of C{x) and D{x) for ^ = 3 and n = 2. The three line segments are the disjoint 
intervals that make up E. The crosses denote the zeroes of C{x) and the dots, the zeroes of D{x). 



is that both sides of (4.24) have the same sign at a; = — 1 + e where e is a small positive constant. This is true 
because close to a; = — 1, Pj{x) and Qj{x) must have opposite sign. Similarly, T„ (^^^^^ and U„ ^ ^a^^^ ) 
have the same sign if R is even and have opposite sign if ^ is odd. The product 11^1^ (a; — /3, ) , for a; = — 1 + e 

is negative if ^ is even and positive if ^ is odd. A^^ is a positive constant, and (1 — x^) is also positive. 
Therefore both sides of (4.24) have the same sign at a; = — 1 + e. 

Therefore if we draw lines connecting the dots and the crosses that represent the tonicity of the polyno- 
mials we will see the pattern shown in figure 9. 

We are interested in where these polynomials intersect, as these points represent the zeros of Pnfi+j{x). 
Starting from the point a; = — 1 and moving along the interval [—1, 1], the following behaviour is observed. 
Before the first occurrence of two adjacent crosses the polynomials intersect between every dot-cross pair. 
After the appearance of the two adjacent crosses, the polynomials intersect between every cross-dot pair. This 
continues until we reach the first occurrence of two adjacent dots and then the intersection of the polynomials 
appears again in between every dot-cross pair. So every appearance of two crosses or two dots swaps the 
position of the intersection points between dot-cross pairs and cross-dot pairs. Since the appearance of two 
dots and two crosses alternate, after every appearance of two crosses there must be a point of intersection 
before the next dot. This corresponds to the statement in the theorem that between every zero of Pj{x) and 



• X — • — X — X- • — X — • — X • •-x — — X — • 

Figure 8: 
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Figure 9: The intersections of the two lines represent qualitatively where the zeroes of PnA+j{x) will lie 
relative to the zeroes of both sides of (4.24) 

the next clement from the set B lies a zero of PnSi+j i^) ■ Similarly, after the appearance of two dots there 
must be an intersection point before the next cross. This is equivalent to saying that between every zero of 
Qj {x) and the next element from set A also lies a zero of Pn^+j (x) ■ □ 

Corollary If Pj{x) has a zero in Ei then so does PnM+j{x)- (Ei := {ai,/3i)) 

Proof: Referring to fig 9 we see that there is always a dot at the beta points. If there is a zero of Pj (x) 
in Ei then the previous theorem tells us that there must be a zero of PnA+j between this zero and the next 
/3 point. Since Ei ends at the next /? point the zero of PnA+j is also contained in £'j. □ 
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Appendix A 

Genus one recurrence coefficients 

Here we show one method of determining the recurrence coefhcients for the genus 1 case explicitly in 
terms of the branch points a and f3. It involves manipulating the theta function expressions given in (1.8) 
and (1.9). We see explicitly how the Jacobian elliptic functions define the form of the coefficients. From [S], 
for g = 1 we have, 



, 1 L ^3("+) , ^!i((2n-3)^+) _ iS',{{2n-l)u+) 
2 A [_ i?3(w+) i?3((2n-3)M+) i?3((2n - l)u+) 



where 



AK 



v/(l-a)(l + /3) 



and. 




with 



2(/3-a) 



(l-a)(l + /?)- 



Using, 



and, 




we write, 



(3 -a 2K l9'o{2Ku+ + K) 9'o{(2n - 3)2Ku+ - K) 0'o{{2n - l)2Ku+ + K) 
2 ^ ~ [ 9o{2Ku+ + K) ^ 6io((2n - 3)2Ku+ - K) ^ 9o{{2n - l)2Ku+ + K) 



The Jacobi function is defined to be 



Z{w) 




and it satisfies the following formula. 



Z{u + v)- Z{u -v)- 2Z{v) 



(M)sn(z;)cn(7;)dn(7;) 



1 — fc^sn^wsn^w 
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If we make the choice, 



u = {2n- 2)2Ki 
V = 2Ku+ + K, 



then we have, 
bn 



2{P-a) 



v/(l-«)(l + /3) 

(3 ~ a 



sn2((2n - 2)2Ku+)ii\i{2Ku+ + K)cn{2Ku+ + K)Ari{2Ku+ + K) 



1 - fc2sii2((2n - 2)2Xu+)sn2(2A'u+ + K) 



Applying some basic identities (fc' — yl— A?), 

cn(ti;) 



and, 



sn(w + K) 



cn(u; + K) = -fc' 
dn(i(; + K) 



I sn(w) 



dn(u;) 
k' 



1 2 

sn + cn w 



dn(w) 



we can write, 



bn = if3- a) 



1 + a 



(/3 - a) - (1 + /3)/sn2((2n - 2)2Ku+) 



smce, 



sn(2i4:u+) 



f3+l 



We have the foUowing addition formula for sn(u) valid for any u and v, 
sn(u)cn(w)dn(u) + sn(w)cn(w)dn(u) 



sn(M + v) 



1 — k^sn^{u)sn^{v) 



We can iterate this to obtain sn((2n — 2)2Ku'^) in terms of k and sn(2i^u+). This gives us in principle, a 
method for evaluating the recurrence coefficients in terms of a and (3. Owing to the nature of the addition 
formula, after just a few iterations the expressions for the coefficients become very large. 



In the a„ are given as. 



yl^9l(2M^ 



Mi'^n + l)u+)M{2n - 3)u+) 
m2n-l)u+) 



71 > 2. 
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Using the following identity, 



1 - k'sn'{u)sn'{v) = 0o(O)' 



with 

u ^ {2n~ l)2Ku+ + K 
V = AKu^ 

we have, 



c (1 - k^sn^{{2n ~ \)2Ku^ + K)%v? {'^Ku^)) n > 2. 



Here c is a constant independent of n, 
/ i?;(0)^9o(2w+ 



c : 



VAz9o(0)??i(2u+) 



To determine c explicitly, we evaluate 02 by comparing coefficients of (1-13) for n = 2 and then compare this 
with a2 given above. In which case we find, 

c=^(2-« + /3)^. 

This value for c was also determined in IHI pg. 4695. Note that in the higher genus cases, the non linear 
difference equations derived in section 2 can be used to evaluate the recurrence coefficients. 
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Appendix B 

g=2 Auxiliary polynomials 

Using the method outlined in section 2 for g = 2 we find, 

G3(x,n) 1 / ^ ^ X 9 

rin-i ^ 



and, 



and. 



hn 



8 

-^iPip2 + a2pi + a2p2 + oi\oi2 + Q!i/3i + ax^i) - 2a„^ x 

+ ^( " (a? + "2 + /3? + + 4(ai + Q2 + /3i + /32) 

+ /3i/32(/3i + /32) + a2/3i(a2 + A) + 0:2/32(0:2 + /32) + Q:ia2(ai + 0:2) 

+ ai/3i(ai + + 0:1/32(0:1 + P2) 

- 2(a2/3i/32 + ai/3i/32 + o:ia2/3i + 0102/32) 

- 16a„(ai + 02 + /3i + /32 - 2(6„ + bn+i))) , n>2 



2x2 - [(ai + aa + /3i + /32) - 26„+i] x 

i (4 + a? + /3? + + /3| - 2(aia2 + ai/3i + a2/3i + ai/32 + a2/32 + /3i/32) 



— 8(a„+i + an) + 46„+i(oi +a2 + Pi+ /32) - 86^_^i) , n > 1. 
^3(2;; 1) = (a; - bi){x - ai){x - 02) 



S2{x; 0) = {x - ai){x - a^)- 
Incidentally, this gives us the 7i,2('T') coefficients for the genus 2 case as, 
7i,2(n) = ^(ai +^2 +/3i +/32) - 26„+i 

T ^8 + 3{al + al+ + /3^) - 2(aia2 + Q:i/3i + a2/3i + ai/32 + 0:2/32 + /3i/32) 

- 16(a„ + a„+i)+46„+i(ai+a2+/3i+/32)-1262_^i ) , n > 1. 
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